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ABSTRACT 

We report on sensitive phase- referenced and gated 1.4-GHz VLBI radio observa- 
tions of the pulsar PSR J0205+6449 in the young pulsar-wind nebula 3C 58, made in 
2007 and 2010. We employed a novel technique where the ^105-m Green Bank tele- 
scope is used simultaneously to obtain single-dish data used to determine the pulsar's 
period as well as to obtain the VLBI data, allowing the VLBI correlation to be gated 
synchronously with the pulse to increase the signal-to-noise. The high timing noise of 
this young pulsar precludes the determination of the proper motion from the pulsar 
timing. We derive the position of the pulsar accurate at the milliarcsecond level, which 
is consistent with a re-determined position from the Chandra X-ray observations. We 
reject the original tentative optical identification of the pulsar by Shearer & Neustroev 
(2008), but rather identify a different optical counterpart on their images, with i?-band 
magnitude ~24. We also determine an accurate proper motion for PSR J0205-I-6449 
of (2.3 ± 0.3) mas yr~-^, corresponding to a projected velocity of only (35 ± 6) km s^^ 
for a distance of 3.2 kpc, at p. a. —38°. This projected velocity is quite low compared 
to the velocity dispersion of known pulsars of ^200 km s~^. Our measured proper 
motion does not suggest any particular kinematic age for the pulsar. 
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1 INTRODUCTION 

The supernova remnant 3C 58 (G130. 7-1-3.1) has a filled cen- 
tre morphology and a flat radio spectrum, and was classified 
on this basis as being pulsar wind nebula (PWN, also ple- 
rion; e.g., Weiler & Panagia 1978) long before a pulsar was 
seen. The pulsar, PSR J0205-(-6449, was subsequently de- 
tected, flrst in the X-ray (Murray et al. 2002) and then in 
the radio (Camilo et al. 2002). It is also one of the approx- 
imately 100 pulsars that shows pulsed gamma-ray emission 
(Abdo et al. 2009). 

PSR J0205+6449 and 3C 58 are at a distance, D, of 
~3.2 kpc (Roberts et al. 1993)^. They have traditionally 
been associated with the historical supernova of 1181 AD 
(SN 1181; Clark & Stephenson 1977; Stephenson & Green 



^ Note that the distance is somewhat uncertain, and somewhat 
larger value is possible. The pulsar dispersion measure is approx- 
imately twice what is expected at the distance of 3.2 kpc (Camilo 
et al. 2002). 



1999), giving them an age of 830 yr and making them one of 
the youngest known supernova remnants and pulsars. How- 
ever, recent work has suggested that 3C 58 is likely consid- 
erably older. In particular, the measured expansion speeds 
of both the synchrotron bubble (Bietenholz et al. 2001; Bi- 
etenholz 2006) and of the thermal filaments (Fesen 1983; Fe- 
sen et al. 1988; van den Bergh 1990; Rudie & Fesen 2007), 
seem to be considerably lower than expected for an age of 
~ 830 yr, suggesting an age of several thousand years, which 
larger age is also suggested several other arguments (see, 
e.g., Bietenholz 2006; Chevalier 2005). 

PSR J0205-I-6449 has a spin frequency of ~15.2 Hz, with 
a derivative of ~ -4.5 x 10"" Hz s"^ (Murray et al. 2002; 
Camilo et al. 2002; Livingstone et al. 2009a), so the char- 
acteristic age of the pulsar is 5400 yr. (We note that if the 
pulsar had an initial spin frequency of ~16.3 Hz, the present 
spin frequency and spindown rate could be reconciled with 
an age of 830 yr; however as mentioned, numerous other 
arguments independent of the characteristic age suggest a 
considerably larger age). It has a high level of timing noise. 
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having exhibited two spin-up ghtches with fractional magni- 
tudes of Aj^/i/ = 3.4 X 10"'' and 3.8 x 10"^ between MJDs 
52276 and 53063 (Livingstone ct al. 2009a). The pulse width 
is ~2.5 ms (at 1.4 GHz, Camilo et al. 2002). Although as 
mentioned, there is some controversy over the age of 3C 58, 
PSR J0205-I-6449 is nonetheless one of the youngest known 
pulsars^ . 

PSR J0205+6449 is a particularly interesting case bo- 
cause of the presence of the easily observable pulsar wind 
nebula, 30 58. Such nebulae, seen for only a handful of pul- 
sars, provide important diagnostics for young pulsars, giving 
insight into the winds which carry away the pulsar spindown 
energies. 

Knowing PSR J0205-|-6449's proper motion is impor- 
tant for determining its exact birthplace as well as the nature 

of its interaction with its surrounding PWN. More generally, 
the origin of the high space velocities of pulsars is an interest- 
ing question, and therefore knowledge of PSR J0205-|-6449's 
proper motion is particularly important because of its young 
age. 



2 PULSAR TIMING WITH GBT: 
OBSERVATIONS AND RESULTS 

We obtained two sessions of pulsar timing observations at 
the ~105-m diameter Green Bank Robert C. Byrd telescope 
at Green Bank (GBT), which were carried out simultane- 
ously with the VLBI sessions. We used a new observing 
technique where the data from the GBT was simultaneously 
used for both pulsar timing and the VLBI observations. The 
VLBI data are described in the next session. 

On 2007 Apr 25, wc observed PSR J0205+6449 with 
the GBT Spigot back-end (Kaplan et al. 2005) for a dura- 
tion of 10.8 hrs. The Spigot measures 3-level autocorrela- 
tions from each of 2 polarisations for 1024 lags and inte- 
grates them for 81.92 (is before summing them and writing 
them to disk. The total bandwidth covered 800 MHz, cen- 
tred at 1650 MHz, approximately 500 MHz of which (from 
~ 1350 — 1850 MHz) were free enough from interference to 
be used for the timing analysis. Wc let the Spigot continue 
to take data while the GBT occasionally moved to and from 
a calibrator source as part of the VLBI scheduling (see be- 
low). For these times, which corresponded to about 35% of 
the total duration, we simply zero-weighted the resulting 
data from the pulsar. A standard pulsar timing analysis us- 
ing TEMPO provided us with a pulsar spin ephemeris which 
we used for pulsar gating of the interferometric data during 
correlation. We measured a barycentric pulsar spin period 
of 0.0657176116(4) s at epoch MJD 54215.3220 (UTO). 

We performed a very similar observation with the newer 
pulsar back-end GUPPI (Ransom et al. 2009) on 2010 Oct 
18, where we observed PSR J0205-F6449 for 8.9 hrs. GUPPI 
recorded the 8- bit, summed polarisations from 2048 channels 
every 192 /is covering 800 MHz of bandwidth (of which the 
highest ^600 MHz was usable). After zero- weighting data 
when the GBT was off of the pulsar's position, a timing 

^ Several pulsars or PWNe, including the Crab Nebula, 
G21.5-0.9 (Bietenholz & Bartcl 2008) and Kcs 75 (Gotthclf et al. 
2000) are thought to be only around 1000 yr old, but the vast ma- 
jority of known pulsars have ages > 10,000 yr. 
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Figure 1. Pulse profiles obtained with GBT at our two observ- 
ing epochs. The profiles have 128 bins, and have been shifted to 
so that the peaks occur at phase = 0.5, and the zero-level in 
each profile is indicated by the dashed horizontal line. We do not 
consider the differences between the two profiles significant. 

analysis dctc>riiiincd the barycentric pulsar spin period to 
be 0.0657388555(6) s at epoch MJD 55487.8302 (UTO). 

For our timing analyses we used a simple two-Gaussian 
model as our template profile where the peak of the pulsar 
flux occurred very near to spin phase 0.5. Wc assumed a 
pulsar position, dispersion measure, and instantaneous spin- 
down rates consistent with those measured by Livingstone 
et al. (2009b). 

We show the average pulse profiles obtained in the two 
observing sessions in Figure 1. As the pulsar is quite weak, 
the signal-to-rioise ratio is relatively low, and consequently 
we do not consider the differences between the two profiles 
significant. The FWHM duty-cycle the pulse is ~6%. 



3 VLBI RADIO OBSERVATIONS 

We obtained two epochs of 1.4-GHz VLBI observations of 

PSR J0205+6449 using the "High Sensitivity Array" , which 
consisted of the NRAO Very Long Baseline Array (VLBA) 
augmented by the GBT (~105 m diameter) and Effelsberg 
(100 m diameter) telescopes. For the first session, on 2007 
April 25 (program code BB241), we obtained 10 hours of 
VLBI data, while for the second session on 2010 October 
18 (program code BB295), we obtained 12 hours. In each 
session we recorded both senses of circular polarisation, and 
we used 2-bit sampling at a bit rate of 512 Mbit s~^, for an 
effective bandwidth per polarisation of 64 MHz. 

The data from 2007 were correlated with NRAO's 
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VLBA processor, while those from 2010 were correlated with 
the DiFX processor (Deller et al. 2011). The analysis was 
carried out with NRAO's Astronomical Imago Processing 
System (AIPS). During correlation, we gated the correlator 
using the pulsar timing extracted from the GBT data, and 
accumulated several sets of VLBI data with different gat- 
ing: firstly, the un-gated data, and secondly, the "on-pulse" 
gated data where the gated bin had a width of 8% of the 
pulsar period (distributed symmetrically about the peak) . In 
addition, for the 2010 data, correlated with the DiFX corre- 
lator, we also accumulated data in 12 different bins evenly 
distributed in pulse phase, with each bin therefore having 
a width of 8.33% of the period. In all data sets we used a 
correlator integration time of 3.0 s. 

We used VCS2 J0209-h6437 (Fomalont et al. 2003), 
0.5° away from PSR J0205+6449, as our primary phase- 
reference source, but wc included some observations of an as- 
trometric check source, ,1VAS J0228-h6721(4C 67.05, ICRF 
J022850.0-I-672103; lERS B0224-H671), which is a source in 
the International Celestial Reference Frame catalogue whose 
position is accurately known (Fey ct al. 2009), and which 
is 3.4° away from PSR J0205+6449. For the phase-centre 
position for the pulsar observations, wc use the position 
of the pulsar given in Slane et al. (2002), 02'^ 05™ 3r!92, 
64° 49' 42'.'8. Mostly we used a cycle time of 3.8 minutes, 
with 2.5 minutes being spent on PSR J0205+6449. The cal- 
ibration of the VLBI data was done using standard proce- 
dures using NRAO's AIPS package. The flux density cali- 
bration was done through measurements of the system tem- 
perature at each telescope, and the antenna amplitude gains 
were subsequently improved through self-calibration of the 
reference sources. 

Some sporadic radio-frequency interference (RFI) was 
seen, particularly in the frequency range 1617 to 1625 MHz. 
Wc clipped out visibility points with anomalously high am- 
plitudes, removing ^$1% of the data. We also deleted any 
data taken when either of the antennas were at elevations 
of < 10° . During our 2007 observing run, the GBT failed to 
observe J0209-I-6437 for three periods of approximately one 
hour. We discarded all the GBT data for all sources from 
these periods. 

The fringe-fitting and phase calibration was done using 
the (un-gated) data from our calibrator sources. We esti- 
mated the ionospheric delay using the AIPS task TECOR 
using lONEX data from the Crustal Dynamics Data Infor- 
mation System archive^ of the Goddard Space Flight Cen- 
ter. The earth orientation parameters used during the cor- 
relation were extrapolated. We corrected to more accurate 
ones subsequently made available from the United States 
Naval Observatory. The resulting calibration was then in- 
terpolated to the times of the PSR J0205-I-6449 as well as 
the J0228-I-6721 observations, and applied to both the gated 
and un-gated visibility data for PSR J0205-I-6449. 



3.1 Astrometric Reference Sources 

Our primary phase-reference source was J0209+6437. We 
take, as the astrometric reference point for J0209-I-6437, the 
peak brightness position, for which we take the coordinates 
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Figure 2. VLBI images of J0209-f6437 on our two observing 
epochs. The contours are drawn at 2, 5, 10, 30, 50, 70, 90 and 
99% of the peak brightness, with the 50% contour being em- 
phasized. The two images have been convolved with the same 
restoring beam of FWHM 4.4 X 3.3 mas at p. a. 8° indicated at 
lower left. The peak brightness was 81 and 59 mjy bm^-' for 2007 
and 2010, respectively. In each panel, the origin is placed at the 
peaJj-brightness point, which we use as the astrometric reference 
point. 
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Figure 3. VLBI images of J0228-t-6721 on our two observing 
epochs. The contours are drawn at 2, 5, 10, 30, 50, 70, 90 and 99% 
of the peak brightness, with the 50% contour being emphasized. 

The two images have been convolved with the same restoring 
beam of FWHM 4.4 X 3.8 mas at p. a. 0° indicated at lower left. 
The peak brightness was 601 and 421 mly bm"^ for 2007 and 
2010, respectively. The visibility data were self-calibrated in phase 
before imaging to best show the intrinsic structure. 



http : / /cddis . gsf c. nasa . gov 
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09™ 35=98806(20), 64° 37' 25'.' 7701(06) from the VCS2 
Survey (Fomalont et al. 2003), where the numbers in paren- 
theses give the uncertainty in the last two digits. Unfortu- 
nately, this source turned out to be somewhat resolved. We 
show the VLBI images of J0209-I-6437 on our two observing 
epochs in Figure 2. The source is elongated in the north- 
south direction, with the 50% contour having an extent of 
approximately 6 mas. 

It is possible that the peak brightness position, and 
therefore our astrometric reference point, varies with time 
due to changes in the source morphology (see, e.g., Bieten- 
holz et al. 2000; Bartcl et al. 2012) 

As mentioned, we also included some observations of an 
astrometric check source, J0228-I-6721. Like our observations 
of PSR J0205-I-6449, we phase-reference our observations of 
J0228+6721 to J0209-I-6437. We give the peak brightness 
positions of J0228-I-6721, as well as the ICRF position in 
Table 1. 

Unfortunately, this source is also extended at the mas 
level in the N-S direction (Lister et al. 2009a; Romney 
et al. 1984; Pcrley ct al. 1980). Moreover, this source is 
known to have a somewhat variable morphology (Lister et al. 
2009b; Bondi et al. 1996; PadrieUi et al. 1986), with both 
changes in the size of the central component, and motions of 
other components being observed with magnitudes of up to 
~ 0.3 mas yr'^ We show images of J0228-I-6721 in Figure 3. 

3.2 PSR J0205-f-6449 

The pulsar was expected to be quite faint: even in the images 
made from gated VLBI data, we expected a signal-to-noise 
ratio of ^10. Furthermore, the pulsar position is not accu- 
rately known a priori. In order not to bias our analysis by 
searching only near an assumed position, we therefore chose 
to image a fairly large region. We imaged PSR J0205-I-6449 
directly from the calibrated multi-channel VLBI data with- 
out any averaging of the visibility data in frequency. Since 
our channel width was 0.5 MHz, the field of view over which 
bandwidth smearing is small is >10", which is adequate for 
our purposes. We made images of 4096 x 4096 pixels, with 
the pixels being 0.8 x 0.8 mas, centred on our initial position 
estimate for PSR J0205-I-6449, which was the position was 
used in the correlation, and was the one given by Slane et al. 
(2002) of 02'' 05™ 37!92, 64° 49' 42'.'8 (J2000). Our images 
therefore span a region ± ~ 2" from this estimated posi- 
tion, which should be larger than any expected uncertainty 
in the X-ray derived position. We used natural weighting 
to obtain the highest signal-to-noise ratio. We evaluated the 
image rms brightness over the central 90% of the image area, 
excluding a narrow region around the edge where the rms is 
typically slightly higher due to gridding artifacts. 

In Table 2 we give the signal-to-noise ratios of the 
brightness extrema in the various images. In the images 
made from un-gated visibility data, the positive and negative 
extrema were of similar magnitude, being ^ 5.5 a where a 
was the image brightness rms. In the images made from the 
on-pulsc gated data, for both epochs, the negative extrema 
were of a similar magnitude. For both epochs however, the 
images made from the on-pulse data showed a single positive 
peak that was ^ 6.5(7, or larger in magnitude than the neg- 
ative extremum by ^ la. The number of independent sky 
positions sampled by our images can by estimated by divid- 
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Figure 5. The inilsc profile of PSR J0205-I-6449: the image 
brightness, given in units of the image rms (<t) at the location 
of PSR J0205-I-6449 in the images made from the data in each of 
12 bins distributed across the pulsar period. The first bin on the 
left is duplicated on the right for clarity. The shaded area shows 
the ±1 o" region. Only in the bin at pulse phase = 0.5 is a signal 
substantially in excess of the noise visible. 

ing the image area by the FWHM area of the fitted beam, 
and was ~650,000. If the pixel brightnesses are Gaussian- 
distributed, the likelihood of obtaining a value in excess of 
±6(T in 650,000 trials is then ~0.1%, although we note that 
the distribution of pixel brightnesses may not be accurately 
Gaussian at these large deviations from 0. 

Furthermore, the two positive extrema seen in the im- 
ages made from the gated data were quite close to one an- 
other on the sky, being within 6 mas of each other on the 
sky. Since the total area we examined was ~ 10^ mas^, the 
odds of random peaks occurring so close together would be 
~ 80000 to 1 (which probability is independent of the actual 
distribution of pixel brightnesses). 

In Figure 4 we show a section of the images made from 
the gated on-pulse and un-gated visibility data sets from 
2010. Although the signal-to-noise is not high, the pulsar 
can be clearly seen in the on-pulse image, but not in the 
corresponding un-gated one. 

Despite the modest signal-to-noise ratio, wo think that 
the detections are firm for several reasons. Firstly, in each of 
the "on-pulse" images, the positive extremum in the image 
(excluding a ~100-pixel strip around the edge where the 
rms is somewhat higher due to numerical artifacts) is ^ 
1 a higher than the corresponding negative one, which is 
not expected to happen by chance. Secondly, no source is 
visible in either of the un-gated images at this location. Any 
background source should be much more strongly detected 
in the images made from the un-gated data since the gated 
visibility data represents only ~8% of the total observations. 

As an additional check, for the 2010 observations, we 
imaged separately the visibility data from the 12 bins spaced 
evenly across the pulsar period. In Figure 5 we show bright- 
ness at the location determined from the on-pulse image as 
a function of the pulse phase. The profile can be compared 
to that obtained from the GBT pulsar-timing observations 
shown in Figure 1 above. 

We note here that position of PSR J0205-I-6449 is dif- 
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Table 1. Position of J0228+6721. 



Observations 


Frequency 
(GHz) 


RA (J2000) 

hr m s 


Decl. (J2000) 

o / // 


Offset ( 
RA 


mas)'' 
decl. 


ICRF (Fey ct al. 2009) 
2007 VLBI obs.t' 
2010 VLBI obs.t- 


8.4 
1.4 
1.4 


02 28 50.05148948 
02 28 50.0516247 
02 28 50.0516134 


67 21 03.0293039 
67 21 03.033364 
67 21 03.035272 


0.78 
0.72 


4.06 
5.97 



''Offset of measured position from the ICRF position. 

''Positions measured from peak-brightness point, and measured relative to J0209+6437, which 
is assumed to be at the VCS2 position (see text). 




MilliARC SEC MilliARC SEC 



Figure 4. VLBI images of PSR .10205+6449 from our 2007 observations. Both images are dirty or undeconvolved. In both panels, the 
grey-scale runs from —3.8 0- (white) to 6.9 u (black), and the FWHM of an elliptical Gaussian fit central part of the beam is shown at 
lower left. The left panel shows the image made from the on-pulse gated visibility data, while the right panel shows the image made from 
the un-gated visibility data. The images are centred on PSR J0205-I-6449, only visible in the left panel, which is at 02*^ 05™ 37!92322, 
64° 49' 41'.'3319. 



Table 2. Position of PSR J0205+6449. 



Observation 

Epoch 


Gating 


Image Extrema 
(in units of image rms) 


RA (J2000) Decl. (J2000) * 
hr m s o ; // 


2007 


Off-pulse 


-5.1 +5.1 


(not detected) 


2007 


On-pulse 


-5.0 +6.8 


02 05 37.92322 (1) 64 49 41.3319 (4) 


2010 


Off-pulse 


-5.5 +5.2 


(not detected) 


2010 


On-pulse 


-5.1 +6.5 


02 05 37.92247 (1) 64 49 41.3343 (4) 



^ The position of the brightness peak of PSR J0205+6449. The digit in parentheses gives the 
statistical uncertainty in the last digit of the corresponding coordinate value; note that the 
total uncertainty is dominated by a systematic component discussed in the text. 
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ferent than that of the point source visible in the X-ray as 
given in Slane et al. (2002) of 02'" 05™ 3r!92, 64° 49' 42'.'8 
by ~ 1'.'5. 

3.3 Proper Motion of PSR J0205-|-6449 

Using the position determinations in Table 2 we obtained 
the displacement of PSR J0205+6449 between our two ob- 
serving epochs to be -4800 ± 570, -1-2800 ± 570 jias in RA 
and DecL, respectively, where the uncertainties are statis- 
tical only. At a distance of 3.2 kpc, the parallax of PSR 
J0205-I-6449 is expected to be 313 fias, smaller than the 
statistical uncertainty, but as the angular displacement due 
to parallax is readily calculable we corrected for it. The 
expected shift between our two particular epochs due to 
parallax is to be 36 and 499 /las in RA and DccL, re- 
spectively. Correcting for the parallax, we then arrive at 
a net displacement of PSR J0205-I-6449 between our two 
epochs of -4860 ± 570, -h2330 ± 570 /ias in RA and DecL, 
respectively (where we have ignored the small additional 
uncertainty incurred due to the uncertainty in the dis- 
tance). This displacement corresponds to a proper motion 
of -1400 ± 160 (Uas yr"^ in RA and 540 ± 160 //as yr"^ in 
DecL, or 1500 ± 160 /las yr~^ at p. a. 159 ± 6°, where the 
proper motion is measured relative to the brightness peak 
of J0209-I-6437, and the uncertainties are statistical only. 

Since our reference source, J0209-I-6437, is somewhat 
resolved, we must carefully assess any possible effect of tem- 
poral changes in J0209+6437 on our proper motion mea- 
surements. In our VLBl observations, we observed a check 
source, J0228-I-6721 in a fashion similar to PSR J0205+6449, 
in particular similarly phase-referenced to J0209-I-6437. The 
positions measured for the check source were given in Ta- 
ble 1, and suggest a proper motion for .10228+6721 of — 1 
and 550 /xas yr~^ in RA and DecL, respectively. The source 
J0228-I-6721 is an ICRF source and Feissel- Vernier (2003) 
finds a proper motion of < 50 /ttas yr~^ at 8.4 GHz over a 
period of ~12 years. However, the lack of secular motions 
at 8.4 GHz, where the core is likely more dominant, does 
not preclude apparent motions at our lower frequency of 

1.4 GHz, where jet components are more likely to dominate 
the emission. Indeed, other authors have reported larger 
proper motions for components within J0228-I-6721: Bondi 
ct al. (1996) find an expansion of 300 /las yr~^ while Lister 
et al. (2009b) find several moving components with proper 
motions of up to 376 /xas yr"'^ at 14 GHz. It is therefore 
possible that our observed proper motion is due to compo- 
nent motions within J0228-f6721. It is, however, also pos- 
sible that our primary reference source, J0209-I-6437, is not 
in fact stable. Unfortunately, both our reference sources are 
elongated in the north-south direction, so that the observed 
north-south relative motion between them could be the re- 
sult of motions along the jet axis in cither source. 

As a conservative estimate of 1 cr uncertainty on the 
proper motion of PSR J0205-I-6449, as phase-referenced to 
J0209-I-6437, we therefore take the apparent proper mo- 
tion of J0228+6721 found above. We then arrive at a fi- 
nal value for the proper motion of PSR J0205-I-6449 of 
-1400 ± 160 /xas yr'^ in RA and 540 ± 575 /ias yr'^ in 
DecL At a distance of 3.2 kpc, this corresponds to a speed 
of 23 ± 6 km s~^. Note, however, that this velocity is with 
respect to the Earth. We calculate a more physical value of 



PSR J0205-|-6449's projected velocity, corrected for Galactic 
rotation and the Sun's motion, in § 5 below. 



4 X RAY ASTROMETRY 

To determine the X-ray position of PSR J0205+6449, we 
investigated data from a deep Chandra observation car- 
ried out between 2003 April 22 and 26 (see Slane et al. 
2004, for details and original results from these observa- 
tions). The data from observation IDs 4383, 4382, and 3832 
were reprocessed and cleaned using standard routines from 
Ciao version 4.4. The merge_all task was used to create 
a merged image from the 317 ks of good exposure time. 
The pulsar is embedded in a bright compact nebula that 
is slightly asymmetric in the cast-west direction. The ccn- 
troid of the point source emission is located at a posi- 
tion of 02'' 05™ 37!93, +M° 49' 41.4'.'. We identified three 
X-ray point sources in the field that have counterparts in 
the 2MASS catalogue. Comparing the centroid positions of 
the X-ray sources with the infrared positions, we find that 
the uncertainty in the Chandra position of the pulsar is 
ctra ~ 0.06 s and Udec ~ 0.12". The systematic error as- 
sociated with contributions from the compact nebula are 
estimated to be slightly larger than these values, but the 
combined uncertainty is ^1 arcsccond cither direction. 

We note that position given above is slightly different 
from that determined in the shorter Chandra observation 
reported by Slane et al. (2002). We have reprocessed those 
earlier data (observation ID 728) as well, and find that the 
position is in excellent agreement with the value given above. 
We conclude that there were small errors in the initial posi- 
tion reconstructions reported by Slane et al. (2002). 



5 DISCUSSION 

Using gated VLBI we imaged PSR J0205-I-6449 in the cen- 
tre of 3C 58. The pulsar was detected. The position of PSR 
J0205-I-6449 we found from the VLBI observations was ~l'.'5 
distant from the originally published one of a compact X- 
ray source seen in Chandra ACIS observations (Slane et al. 
2002). A re-examination of the Chandra data, however, re- 
sulted in an improved position of the X-ray source which is 
well within the uncertainties of that measured with VLBI. 

Shearer & Neustroev (2008) report of deep optical ob- 
servations of the centre of 3C 58 using the 4.2 m William 
Herschel Telescope in La Palma, in which they detect three 
unresolved sources with 7?-band magnitudes of ^24, in ad- 
dition to the more diffuse synchrotron emission from the 
pulsar wind nebula. Shearer & Neustroev (2008) tentatively 
identified their source "ol" as PSR J0205+6449 on the basis 
of its coincidence with the published position of the com- 
pact X-ray source. However, based on our VLBI determina- 
tion and the re-determined Chandra HRC position, it is in 
fact their object "o2", which is almost certainly the optical 
counterpart of PSR J0205+6449. The position of o2 is RA = 
02'' 05™ 37^93 and DecL = 64° 49' 4l'.'4, with an uncertainty 
of ^O'.'l, consistent to within 0'.'08 with the pulsar position 
determined from VLBI. Shearer & Neustroev (2008) give the 
magnitudes of o2 as 24.15±0.07 in R, > 24.3 in V and < 25.6 
in B, consistent with the optical magnitudes estimated for 
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Figure 6. A 1.4-GHz VLA radio image of 3C 58, reproduced from Bieteniiolz (2006). The "x"sign shows the present location of PSR 
J0205+6449, while the two "+" signs show the extrapolated position at two possible dates of the supernova explosion which gave rise to 
3C 58: at 1181 AD and at 5000 BC, with the one at 1181 AD being the one overlapping with the present position (x). The "o" shows 
the position of the centre of the circular region of softer X-ray emission, interpreted as thermal emission from the supernova shell, from 
Gotthelf et al. (2007). The grey-scale is labelled in mjy bm~^, and the FWHM resolution was 1'.'4. We show a detail of the central region 
in Fig. 7 below. 



PSR J0205-I-6449 based on its spindown luminosity, distance 
and the expected Galactic extinction. 

In summary, we can regard the detection of PSR 
J0205+6449 in the radio, optical and X-ray bands as firm, 
with the most accurate position being that determined from 
the VLBI observations and given in Table 2 above. 

In addition to determining the position, we also deter- 
mined the proper motion of PSR J0205+6449 of -1400 ± 
160 /ias yr~^ in RA and 670 ±575 /^as yr~^ in Decl. In order 
to correct this measured value for Galactic rotation and the 
Earth's motion with respect to the local standard of rest, we 
take the Galactic constants from Schonrich (2012), namely 
a flat Galactic rotation curve with v = 238 km s~^ and 
a distance to the Galactic centre of 8.27 kpc. We take also 
from Schonrich (2012) the Solar motion of 250 km s~^ in the 
Galactic plane, solar radial velocity of 14 km s^^ towards the 
Galactic centre, and a motion of 6.1 km s~^ perpendicular to 
the Galactic plane. We take the distance of PSR J0205-I-6449 
(from Earth) to be 3.2 kpc. We can then calculate the pecu- 
liar motion of PSR J0205-I-6449 with respect to the standard 
of rest at its location as being* 2.3 ± 0.3 mas yr~^, corre- 



sponding to 35 ± 6 km s~^ at p. a. —38° (p. a. in equatorial 
rather than Galactic frame). 

This speed is relatively small: Galactic pulsars have a 
2-D velocity dispersion of 200 ~ 300 km s~^ (e.g., Hobbs 
et al. 2005; Lyne & Graham-Smith 1990). Even taking the 
uncertainty introduced by possible motion of our reference 
source, the 3cr upper limit on the proper motion of PSR 
J0205+6449 is only 53 km s~\ 

Can PSR J0205+6449's low tangential velocity be rec- 
onciled with the velocity distribution of other Galactic pul- 
sars? Faucher-Giguere & Kaspi (2006) found that the dis- 
tribution of space velocities of Galactic pulsars at birth was 
consistent with one where each of the three orthogonal com- 
ponents of the space velocity was distributed exponentially, 
with the mean of the absolute value of each component being 
180^30 km s~^, and the resulting mean three-dimensional 
speed being 380tg|5 km s"\ Since PSR J0205+6449 is quite 
young (we discuss the age in more detail below), its ve- 
locity is probably close to its birth velocity. As the multi- 
dimensional exponential distribution of Faucher-Giguere & 



We noted above that the proper motion uncertainty is rather 



larger in Decl. than in RA. Here we take the geometric mean of 
the two values. 



© 2011 RAS, MNRAS 000, 1-10 



8 Bietenholz et al. 



[h] 



1.0 1.2 1.4 1.6 1.8 




02 05 40 39 38 37 36 35 34 
Right ascension (J2000) 



Figure 7. Detail of the 1.4 GHz VLA radio image of 3C 58 from 
Figure 6 above. The symbols have the same meanings as in that 
figure. The restoring beam FWHM of 1'.'4 is indicated by the 
circle in the lower left corner. 

Kaspi (2006) is difficult to evaluate numerically, we per- 
formed a Monte-Carlo simulation with n = 10000 trials, 
to find that the chance of finding a pulsar with a tangential 
velocity as small as the 35 ± 6 km s~^ measured for PSR 
J0205-|~6449 from such a distribution is ~2.7%. 

Alternatively, Hobbs et al. (2005) collected proper mo- 
tion measurements for 140 pulsars^. Of their sample, 8, or 
5.7% had tangential velocities 35 km s~^. Finally, some 
authors have found the distribution of pulsar velocities to 
be bi-modal. For example, Brisken et al. (2003) find that 
20% of pulsars form a low-velocity component with a one- 
dimensional velocity dispersion, (j„id of 99 km s~^, while 
the remainder have a o"„id ~ 294 km s"\ In this case, the 
probability of a random pulsar having a tangential velocity 
as low as 35 km s~^ is 1.8%. 

In conclusion, the small measured proper motion makes 
PSR J0205-(-6449 somewhat unusual in being amongst the 
slowest few percent of young pulsars regardless of whether 
a bimodal or a unimodal distribution of pulsar velocities is 
considered. 

Our estimate of the tangential velocity depends on the 
distance to PSR J0205-)-6449 and 3C 58, which, as men- 
tioned earlier, is somewhat uncertain. The value of 3.2 kpc of 
Roberts et al. (1993) was determined kinematically from HI 
absorption. Kothes (2010) argues for a distance of ~2 kpc. 
Adopting a smaller distance would reduce the tangential ve- 
locity estimate by the corresponding factor, and thus make 
PSR J0205-|-6449's low speed even more unlikely. The low 
measured angular speed therefore argues against a distance 
much lower than 3.2 kpc. 

On the other hand, PSR J0205-|-6449's dispersion mea- 
sure is approximately twice that expected for a distance 
of 3.2 kpc (Camilo et al. 2002), suggesting perhaps that 
the true distance is somewhat larger, which would imply a 

^ In fact, Hobbs et al. (2005) collected proper motion measure- 
ments of 175 pulsars, but we exclude the 35 recycled pulsars from 
Hobbs' sample as they are much older and have a notably lower 
velocity dispersion than the remainder. 



tangential velocity larger by a factor of (D/3.2 kpc). This 
question could be resolved with a measurement of PSR 
J0205-|-6449's trigonometric parallax with 10% accuracy. 
Such a measurement is feasible if an in-beam calibrator can 
be found (see e.g., Chatterjee et al. 2009), and should be 
undertaken. 

Can we extrapolate back from PSR J0205-|-6449's 
present position to determine its position at the time of the 
supernova explosion? To do so requires knowing the age of 
PSR J0205+6449. As noted in the introduction, although it 
has traditionally been associated with a supernova in 1181 
AD, making it ~830 yr old, it is probably older, with a likely 
age of ~7000 yr (Bietenholz 2006; Bietenholz et al. 2001), 
suggesting an epoch of ~5000 BC for the supernova event. 
In Figure 6 we show a radio image of 3C 58, and indicate 
the present, and extrapolated epoch 1181 AD and 5000 BC 
positions of PSR J0205-H6449. 

From XMM-Newton X-ray observations of 3C 58, Got- 
thelf et al. (2007; see also Bocchino et al. 2001) found an 
approximately circular region with a softer X-ray spectrum, 
which they interpreted as a thermal X-ray emission from 
the supernova shell. They noted, however, that the cen- 
tre of this shell was displaced from the pulsar position at 
Q2^ 05™ 33=97, 63° 49' 50'.'0 (J2000.0). We plot this position 
also in Figure 6. Neither the present distribution of syn- 
chrotron emission or thermal filaments suggests any par- 
ticular explosion centre, and although Fesen et al. (2008) 
determined proper motions for various optical features, the 
precision is too low to accurately identify an expansion cen- 
tre. 

Our measured proper motion does not place PSR 
J0205-I-6449 near the centre of the region of softer X-ray 
emission for either possible explosion epoch. For an explo- 
sion epoch of 5000 BC, however, the extrapolated position 
of PSR J0205-I-6449 is closer to the geometrical centre of 
3C 58. 

If the area of softer X-ray emission identified by Got- 
thelf et al. (2007) were thermal X-ray emission associated 
with the forward shock of the supernova, then one would 
probably expect its centre to be near the location of the 
explosion and of PSR J0205-|~6449's birth. This does not 
seem to be the case. Furthermore, the diameter of the re- 
gion of softer X-ray emission is only ~5.6 pc, which is no- 
tably smaller than the PWN, which has an E-W extent of 
~8.5 pc. The PWN, however, is expected to still be confined 
by the ejecta and thus be inside the forward shock, which 
suggests that the forward shock is considerably larger than 
the region of softer X-ray emission. Therefore, both because 
the region of softer X-ray emission is not centred on the lo- 
cation of the explosion, and because it is small compared 
to the PWN, we think it is unlikely that the softer X-ray 
emission is associated with the supernova forward shock. 



6 SUMMARY AND CONCLUSION 

1. We obtained VLBI observations of PSR J0205-)-6449, the 
pulsar in 3C 58. We employed a novel technique to obtain 
VLBI observations of faint pulsars, where we used the GBT 
simultaneously for pulsar timing observations as well as an 
element of the VLBI array. The derived pulsar timing infor- 
mation was then used to gate the VLBI correlator increasing 



© 2011 RAS, MNRAS 000, 1-10 



The Proper Motion of PSR J0205+6449 in 3C 58 9 



the signal-to-noise of the pulsar VLBI. This technique can 
be used to advantage for other young pulsars which have 
high timing noise. 

2. We determined an accurate position for PSR J0205-I-6449 
of 02*" 05™ 37!92 64° 49' 41'.'3, which is l'.'5 different than 
the previously accepted one, which was based on an X-ray 
image. Re-examination of the X-ray data, however, reveal an 
error in the original reported position; the newly-determined 
X-ray position reported here is consistent with the VLBI 
position. Furthermore, this position is coincident with an 
optical source identified by Shearer & Neustroev (2008). 

3. We determined the proper motion of PSR J0205+6449. 
After correction for Galactic rotation, we found a proper 
motion of 2.3 ± 0.3 mas yr~^, corresponding to tangential 
velocity of (35 ± 6)(Z)/3.2 kpc) km s"^ at p.a. -38°. This 
low speed puts PSR J0205-I-6449 amongst the slowest few 
percent of young pulsars. 

4. We estimated PSR J0205+6449's position at birth. If it 
is the remnant of a supernova in 1181 AD, then its position 
at birth is only l'.'9 different than at present. If as seems 
more likely, the age of PSR J0205+6449 and 3C 58 is sev- 
eral thousand years, then its position at birth was near the 
midpoint of the presently visible nebula, although somewhat 
displaced from the present radio brightness centre which is 
close to the present location of PSR J0205-I-6449. 
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